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ABSTRACT 

I conduct an analytical study of the interaction of jets, or a coUimated fast 
wind (CFW), with a previously blown asymptotic giant branch (AGB) slow wind. 
Such jets (or CFWs) are supposedly formed when a compact companion, a main 
sequence star or a white dwarf, accretes mass from the AGB star, forms an accre- 
tion disk, and blows two jets. This type of flow, which is thought to shape bipo- 
lar planetary nebulae (PNe), requires 3-dimensional gas dynamical simulations, 
which are limited in the parameter space they can cover. By imposing several 
simplifying assumptions, I derive simple expressions which reproduce some basic 
properties of lobes in bipolar PNe, and which can be used to guide future numer- 
ical simulations. I quantitatively apply the results to two proto-PNe. I show that 
the jet interaction with the slow wind can form lobes which are narrow close to, 
and far away from, the central binary system, and which are wider somewhere 
in between. Jets that are recollimated and have constant cross section can form 
cylindrical lobes with constant diameter, as observed in several bipolar PNe. 
Close to their source, jets blown by main sequence companions are radiative; 
only further out they become adiabatic, i.e., they form high-temperature low- 
density bubbles that inflate the lobes. This implies that radiative cooling must 
be incorporated in numerical codes intended to study the formation of lobes in 
PNe. 

Subject headings: circumstellar matter — planetary nebulae: general — stars: 
AGB and post-AGB — stars: mass loss — ISM: jets 
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1. INTRODUCTION 

Planetary nebulae (PNe) and proto-PNe possess a rich spectrum of different structures, 
with spherical and bipolar PNe at the two extremes of that spectrum. Bipolar PNe are 
defined as extreme aspherical PNe having two lobes with an 'equatorial' waist between them 
(Schwarz, Corradi & StangheUini 1992). Several models have been suggested to explain the 
formation of lobes. One of the popular models a decade ago was the interacting wind model, 
where a fast wind blown by the central star of the PN is coUimated by a previously ejected 
AGB stellar dense equatorial wind (e.g., Balick 1987; Soker & Livio 1989; Frank & Mellema 
1994; Mellema & Frank 1995; Mellema 1995). When the equatorial to polar density ratio is 
very high, a bipolar nebula was supposed to be formed. However, the numerical simulations 
did not form bipolar PNe with very narrow waists, but rather formed elliptical PNe, or 
bipolar PNe with wide waists. Another problem of the interacting wind model is the high 
momentum and kinetic energy observed in several bipolar PNe and proto-PNe (Bujarrabal 
et al. 2001; see Bahck 2000 for some other problems). Therefore, the simple interacting wind 
model cannot explain many of the observed bipolar PNe. It seems that in order to form many 
bipolar PNe, e.g., those with very narrow waists, a coUimated fast wind (CFW) is required 
(Morris 1987; Soker & Rappaport 2000, hereafter SROO). Such preliminary simulations, but 
where the CFW was blown by the primary star, were performed by Frank, Ryu, & Davidson 
(1998), who indeed got very narrow waists in their simulated nebulae. Other types of models 
are based on magnetic fields playing a dynamical role in the mass loss process from AGB 
stars (e.g., Garcia-Segura 1997; Garcia-Segura & Lopez 2000; Blackman et al. 2001). There 
are some severe problems with dynamical magnetic models (Soker & Zoabi 2002). The main 
problem is that the AGB has to be spun-up by a stellar companion in order to possess the 
required activity; hence we are in the regime of binary models. 

As reviewed recently by SROO, the most promising mechanism for the formation of 
bipolar PNe with narrow waists is based on an accreting compact companion to the mass 
losing AGB star. In most cases, the companion is a white dwarf or a main sequence star 
(SROO). When a compact binary companion accretes matter with sufficiently high specific 
angular momentum, an accretion disk is formed. When the accretion rate is above some 
threshold (SROO), two jets are expected to be blown by the companion. If the jets are not 
well coUimated, the outfiow from the companion is termed CFW. In addition, the companion 
may lead to a concentrated equatorial fiow (Mastrodemos & Morris 1999), further increasing 
the equatorial to polar density ratio. The formation of bubbles (or lobes) by jets have been 
studied in a variety of astrophysical fiows, e.g, young stellar objects (e.g., Masson & Chernin 
1993), clusters of galaxies (e.g., Reynolds, Heinz & Begelman 2001), and gamma ray burst 
environments (Ayal & Piran 2001). However, the condition of jets blown into AGB winds 
are different from those in these cases. 
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To study the interaction of the two winds blown simultaneously, the slow wind by the 
AGB star and the CFW (or two jets) by the accreting compact companion (see figs. 1 and 2 
in SROO), a 3D gasdynamical numerical code is required. Such numerical simulations require 
heavy computer resources, hence can cover only a small number of cases. Therefore, it will 
be helpful to derive simple expressions than can both describe the general structure of lobes 
obtained by different jets, and guide future numerical simulations. This paper is aimed at 
reaching these objectives, by exploring part of the parameter space of jets interacting with 
spherical AGB winds, and under simplifying assumptions. To manifest analytical solutions, 
I consider two extreme cases, which are described in §2, and are studied in §3 and §4. A 
short summary is in §5. 



2. JET PROPAGATION 

In the present analytical study, I treat the flow at large distances from the binary system, 
neglecting the deflection of the jets (or CFW; see SROO). I assume an axisymmetric flow 
in which a jet expands along the symmetry axis and into a previously ejected spherically 
symmetric slow AGB wind. The slow wind density as a function of distance from the binary 
system is pg = Ms/{47rVsr'^), where Vg is the slow wind speed, and Mg is the mass loss 
rate, both are constant, and where mass loss rate is defined positively. The gas inside the 
undisturbed jets, one at each side of the equatorial plane, expands at a constant speed Vj, 
the mass loss rate into each jet is Mj, and each jet expands into a solid angle of flj = 47r/?. 
The density inside an undisturbed radially expanding jet, i.e., its cross section increases as 
r^, is pj = Mj/{A7r/3vjr'^). The velocity of the jet head, Vh, is determined by equating the 
slow wind pressure on the jet head with that of the jet material. Since thermal pressure can 
be neglected in the pre-shocked media, the expression is Pj{vj — v^f — Psi^h ~ Vg)^. In the 
present analytical study I treat two extreme cases: slow-propagating and fast-propagating 
jets. In the slow-propagating jet case the jet head proceeds at a speed Vh <^ Vj, while in the 
fast-propagating jet case the jet head does not slow much, Vg -^Vh Vj. 

The velocity of a slow-propagating jet to be treated in the next section is given by 

1/2- 



Vh 



1 + 



(3MgVg 



for pg > Pj. (1) 



I therefore assume that MjVj < few x lOPMgVg, so that Vh ^ 2 — lOvg. Since I neglect 
the bending of the jet by the AGB slow wind close to the binary system, the momentum 
flux of the jet must be (SROO) MjVj > MgVgtana, where a is the half opening angle of 
the jet (SROO). For tana < 1, tana ~ 2(/3)V2. For a = 10°, for example, /3 = 0.0076, 
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2/3^/^ — 0.174, and tana = 0.176. Over all, the treatment of the slow-propagating jet is 
applicable to the case {MjVj)/{MsVs) ~ 0.1 — 0.5, for a half opening angle of a ~ 5 — 30°. 
However, since, as we'll see later, the derived expressions are not so sensitive to the different 
variables, the results will be applicable to a much larger parameter space. 

Prom the expressions developed above, we can find the density of the jet shocked mate- 
rial, and from that, the radiative cooling time, tcooi = (5/2)nA';T/(nenpA), where n, n^, and 
Up, are the total number, electron, and proton, post-shock densities, respectively, taken to 
be four times the pre-shock densities, and for the considered temperature range, which is 
determined by the speed of the jet material, A ~ 10^^^ erg cm^ s~^. The cooling time of the 
shocked jet material is obtained by substituting typical values 

This time scale is much longer than the flow time along the jet diameter. However, the 
cooling time of the post shocked material should be treated differently and compared with 
another flow time scale. This is done for the two flow cases in the next two sections. 



3. SLOW-PROPAGATING JETS 

3.1. Jet Interaction with the Slow Wind 

The shocked jet material will expand and form a hot low density bubble that will 
accelerate the dense slow wind's material around it. The results of Castor, McCray & 
Weaver (1975) for a stellar wind bubble in a dense interstellar medium can be used. Using 
their equation (6) for the radius of the expanding bubble as a function of the expansion time 
of the bubble, th, gives 

o\ 1/5 / . \ 1/5 

0.5M,^;n _ / V, N2/5/ M, 



Rh = 0.76\ ^ f^/^' = 1.0 X lO^M -] = — ^ — -r. 

y Ps J V400 kms-V \^10-^MoyrV5 

X ( ^ ( ^ (^—] cm (3) 

Vl5kms-V llO-^Moyr-M VlO^e cm/ \im ) ' ^' 



where the local slow wind density at r was taken in the second equality. The expansion 
velocity of the bubble's surface is 



X 1/5 

2/5 / Mi \ r V, \V5 



0.6i?fe/t = 20 (^4QQ]^^g_ij I io-7Mq yrV5 I 1 15 km s-V 
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J VlOOvr/ ^ ^ 



lO-SMoyr-i/ VlQi^ cm/ VlOOyr 



The expansion time of the hot bubble, tf,, should be compared with the cooling time of 
the hot gas inside the bubble. The density inside the bubble is given by ph — Mjth/ (47ri?^/3), 
from which the cooling time is derived 

-1 



-1 



t.o.,(bubble) c 30 [ ^^^ j (^^^TO^^ J 1, 100^ J ^"^'^ 

The condition for the hot gas (the shocked jet material) to inflate a bubble is that the cooling 
time is longer than the expansion time, i.e., icooi(bubble) > ti,. Substituting for icooi (bubble) 
from equation (5), and eliminating i?f, by equation (3), gives the radius outward to which 
the shocked jet material will not have time to cool before inflating a large bubble (which will 
evolve to a nebular lobe) 

> 4.4 X 10- ( ) { \ " ( ) 

V400 km s-V I IQ-^M,., vr-^ / Vl5 km s-V 

cm. (6) 




100 yr 

The scaling in the last equation shows that the adiabatic phase, where the cooling 
time of the shocked material is very long, starts at r ~ 10^^ cm. However, the parameters 
change from one system to another. For example, in the proto-bipolar PN OH 231.8+4.2, 
the mass loss rate is very high, 2 x 10~^Mq yr~^ (Alcolea et al. 2002 and more references 
therein). If we take the mass loss rate in both the AGB wind and the jet to be 10 times 
higher than the values used in the scaling of equation (6), we find the adiabatic phase to 
start at Ts ~ 5 X 10^® cm from the central source. If the companion accretes via a Roche 
lobe overfiow, then the mass loss rate into the jets can be much higher. In OH 231.8+4.2, 
the mass loss rate into each of the two jets could have been ~ IO^^Mq yr~^ for a short 
time during a Roche lobe overflow (Soker 2002). If we assume that the opening angles of 
the two jets were large, so the jets' density is lower than the slow wind density, we can use 
equation (6). Substituting M, = 10'^ Mq yr-\ Mj = lO'^M© yr-\ and Vj = 500 km s^S 
we find from equation (6) Vs — 5 x 10^^ cm. This is a substantial fraction of the lobe size 
(Kastner & Weintraub 1995). Therefore, during a large fraction of the evolution (since the 
jet phase was short; Soker 2002) the fiow was momentum conserving (radiative), i.e., the 
shocked material was cooling on a relatively short time. Indeed, Bujarrabal et al. (2002), 
ruled out the adiabatic case for OH 231.8+4.2, mainly on the ground of the shape of the 
lobes, which are not infiated as is expected in the adiabatic fiow. 
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3.2. The Bubble Shape 

From equation (3) it turns out that the radius of a spherical bubble increases with r 
much slower than r does. For example, at a distance from the center 32 times that was is used 
in equation (3) the bubble radius will be only 4 times larger, keeping the time and the other 
parameters unchanged. However, the bubble will not be spherical since the jet propagates 
outward, and a more accurate treatment is required. I assume that the bubble expands 
perpendicular to the symmetry axis, along the z direction, i.e., parallel to the equatorial 
plane. I also take the density in planes parallel to the equatorial plane to be the same as 
the density where the plane intersects the symmetry axis. The distance of the plane from 
the equatorial plane is the coordinate x. The relevant quantity is the energy deposited per 
unit length along the jets propagation, which is cq = 0.5Mjf|/f/i. I use the relation between 
pressure and energy and the momentum equations as in Castor et al. (1975), and neglect 
the thermal energy of the swept slow wind material. The last assumption is adequate since 
I neglect the slow wind radial expansion and the dependence of density on the distance from 
the symmetry axis, while the assumptions allow a simple expression for the energy equation. 
The energy equation reads e — cq — Ck, where e is the thermal energy per unit length (along 
the symmetry axis) in the bubble, and is the kinetic energy per unit length of the swept up 
slow wind gas. Neglecting the density variation perpendicular to the symmetry axis and the 
radial expansion of the slow wind gives the distance of the bubble surface from the symmetry 
axis 



Zb{x) 



16 ep 
57rps(x) 



1/4 



^1/2 



eo 



Ps{x) 



1/4 



(7) 



where x and z are the coordinates parallel and perpendicular to the symmetry axis, respec- 
tively. Substituting in the typical values used here, with Vh by equation (1), gives 



— ~ 1.3 

X 



400 km s" 



1/2 



Mi 



Vh ^ 



yr-1 



X 



yr-1 



tbjx) 
100 yr 



1/2 



X 



-1/4 



-1/2 



10^6 cm/ 



(8) 



For the same parameters, the velocity is 



_ dzb Zb 
Vb = —r- — 0.5— ~ 20 
dt t 



h{x) 
100 yr 



-1/2 



{ttJk- ) km s ^ 

VlO^s cm/ 



(9) 



Because of the assumption that density does not depend on the z coordinate, the last two 
equations are less accurate for Zb > r, i.e., close to the central binary system. 
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It should be noted that in the last two equations the time tb{x) is counted from the 
moment the jet reaches the point x, hence ti is longer for lower values of x, i.e., closer to the 
center of the nebula. If the nebula is observed at a time tobs after the jets are launched, and 
the jets proceed at velocity Vh, then th{x) = tohs ~ ^hh- Substituting this in equation (8), 
keeping the other variables the same, gives 



Zb ~ 1.3 X 10^^ 



(10) 



100 yr J V 1016 cm 
The maximum width of the lobe, ^bmax, occurs at 

. = 1.6 X 10- cm (^) {^^^) cm. (11) 

assuming the jet life time was > tobs/2. For this time of observation, Vh — 100 km s^^, and 
with the other values in equation (8), the maximum width is z/x ~ 0.7. Due to the decrease 
of the density with distance from the symmetry axis and the radial expansion of the slow 
wind, the lobe will actually be somewhat wider. These values for the location of the widest 
lobal point and its width there are typical for many young bipolar PNe. The derivation 
above shows that the bubble, or lobe, is narrow close to the central system, then becomes 
wider, and then narrows again. If the density of the material on the outer portion is too low 
to be observed, the observed lobe surface will only widen with distance from the equatorial 
plane. 

There are several strong assumptions in the derivation above; in particular, (i) neglecting 
the variation of the slow wind density perpendicular to the symmetry axis; {ii) treating 
the slow wind as a static medium rather than a radially expanding medium; and {Hi) the 
assumption that the bubble surface only expands perpendicular to the symmetry axis (the z 
direction). The last assumption is reasonable at the intermediate values of x, but not close to 
the equatorial plane (small values of x) . The expansion there is more appropriately described 
by equation (3) and (4). Also, from equation (6) it is clear that the center of each bubble 
will be at a large distance from the equator. Hence, the expanding bubble surface facing the 
equatorial plane will push material toward the equatorial plane, substantially increasing the 
density there. This supports the claim of SROO that jets (or a CFW) can compress the gas 
near the equatorial plane, contrary to some claims that the dense equatorial gas collimates 
the jets. If the bubble surface reaches the equatorial plane, it will collide with the bubble 
on the other side of the equatorial plane, and the lobe sides will be observed to emerge from 
the equatorial plane at a large distance from the center. Such a structure is observed, for 
example, in the PN He 2-104 (Sahai & Trauger 1998). The inflation of the bubble away 
from the central system implies that the bubble boundary may not reach the equatorial 
plane. This effect is stronger if we consider the outward motion of the slow wind, which 
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was ignored in the derivation above. However, the distance of the bubble surface from the 
equatorial plane will be very small, since the bubble is starting to be inflated by the jet (or 
CFW) at r ~ 400 AU (eq. 6). To observe a bubble not touching the equatorial plane, the 
nebula must be observed almost edge on. These two constraints imply that only a very small 
number of such bipolar PNe should be observed. Taking into account its inclination, it seem 
that the lobes of the bipolar-PN Hb 12 do not cross the equatorial plane (Sahai & Trauger 
1998). At larger values of x the jets may cease to exist, and a different treatment is required. 
Despite these approximations, the treatment above gives the basic shape of lobes in some 
bipolar-PNe, and demonstrates the capability of the binary-blown jets to form bipolar lobes. 

The scaling of the jet velocity, 400 km s^^, is typical for jets blown by main sequence 
stars, which have similar escape velocities (Livio 2000). Jets blown by white dwarfs will have 
much higher velocities, of the order of the escape velocity from white dwarfs (~ 5000 km s~^). 
Such velocities are directly observed in some PNe and proto-PNe, e.g., a wind velocity of 
2300 km s"-*^ was reported by Sanchez Contreras & Sahai (2001) in the proto-PNe He 3-1475. 
In that case, we find from equation (6) that the shocked jet material will not cool, and the 
bubble will be inflated as soon as the jet encounters the slow wind gas. In these cases, the 
lobes that are formed will have a wide opening very close to the equatorial plane, e.g., as in 
the symbiotic Mira He 2-104 (the Southern Crab; Corradi et al. 2001). 



4. FAST-PROPAGATING JETS 

4.1. Jet Interaction with the Slow Wind 

When the density of the gas in the jet is much larger than that of the slow wind it 
will not slow much, and equation (1) is not applicable. Instead, the jet head proceeds at a 
velocity given by 

1/2" 

for Pj > Ps. (12) 



1 - 



ps 

Most of the thermal energy released is that of the shocked slow wind material, rather than 
that of the shocked jet material as in the slow-propagating jet case. Thermal energy injected 
to the shocked slow wind material per unit time and unit length are Eq — 0-5AjPsVj, and 

eo = Eo/vh ^ 0.5AjPsv], (13) 
respectively, where Aj is the cross section of the jet. 

As in the slow-propagating jet, bubble inflation starts when the shocked material has no 
time to cool. To flnd this distance Tg, the calculations that led to equation (6) are repeated 
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with the following changes: (i) the energy injection rate O.bMjVj in equation (1) is replaced 

by Eq given above with Aj = 47r/?r^, and (ii) the shocked material now is the slow wind gas. 
Hence the term for the total shocked mass Mjtb in the derivation of equation (5) is replaced 
by the mass of the shocked slow wind during time t^: jSMgiiiVj/vs, where, as before, is 
the expansion time of the bubble. The condition that the cooling time be longer than the 
expansion time 4, becomes a condition on the location of the jet head 

r > 6 1 X 10^^ ( Y'^' ( r'^' ( ' 

'~ ■ V400 kms-V Vl5kms-V yiQ-^ Mq yr'^ 

X f ^V'V (14) 

VlOOyry/ \lO-y ^ ' 

This condition is very similar to the one for slow-propagating jets (eq. 6) . As with slow- 
propagating jets, jets blown by main sequence stars, for which Vj ~ 400 km s~^, will start to 
inflate a bubble only after propagating a significant distance from their origin; this will push 
slow wind material toward the equatorial plane. Jets blown by white dwarfs companions, 
for which Vj » 1000 km s~^, on the other hand, will start to inflate a bubble very close to 
their origin at r ~ 10 AU. They will also push material toward the equatorial plane, but in 
a complicated flow which required 3D numerical simulations to be treated correctly (SROO). 



4.2. The Bubble Shape 

I consider two cases which bound the two extremes of a jet's possible cross section. In 
the first case, the jet expands radially so that Aj — 47r/3r^, at all times. In the second case, 
the jet expands radially at early times, but it is then recoUimated by the pressure of the 
surrounding shocked gas or by an internal magnetic field (e.g., Kossl, Miiller, & Hillebrandt 
1990), such that its cross section stays constant at Aj = 7ra|, where aj is the radius of the 
jet. 1 do not study these processes here, and simply assume that after opening up, the jet 
is rcconfined near the place where bubble inflation starts, ^ 10^^ cm by equation (14), 
and I scale the jet radius by aj = O.lr^ ~ 10^^ cm. Jets propagating with their cross section 
increasing at a rate much slower than Aj oc r^, and even with constant cross section along 
a fraction of their path, are observed in different objects. Examples are the radio jet in the 
radio galaxy 3C 111 (Linfield & Perley 1984) and the molecular outfiows in the young stellar 
objects HH 288 (Gueth, Schilke, & McCaughrean 2001) and HH 34 (Reipurth et al. 2000). 
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4-2.1. Freely expanding jets 



Substituting the relevant cross section, Aj = An/Sr'^, in the rate of thermal energy re- 
leased per unit length (eq. 13), and then in equation (7), gives the radius of the axisymmetric 
bubble as function of the distance, x, from the equatorial plane 



Zb 

x 



~ 1.0 



400 km s" 



^1/2 


tb{x) 




100 yr_ 



1/2 



10^^ cm 



-1/2 



1010-3 



1/4 



(15) 



We find that the basic shape of the lobe is the same as in the case of a slow-propagating jet, 
hence, the analysis from equation (8) to the end of the previous section is applicable to the 
present case as well. 



4-2.2. Reconfined jets 

Substituting the constant cross section in the thermal energy released per unit length 
(eq. 13), and then in equation (7), gives the radius of the axisymmetric bubble formed by a 
reconfined jet 

1/2 

cm. (16) 

The shape here is different from that in previous cases; rather than widening, the cylindrical 
lobe has a constant radius. Since the jet propagates outward, the expansion time at location 
X from the equatorial plane is t^bs ~ ^/'^j-i where ^obs is the observing time measured from 
the birth of the jet . Since in the present case the jet propagate much faster than the slow- 
propagating case, this delay is less significant than the case analyzed in equations (lO)-(ll). 
In any case, for the same parameters as in the last equation, the bubble radius will be 

~ 4 X 10^6 ( ( \ 

^,_4.UxiU l^4QQ]^^g_i; 1,1015 cmJ 

1/2 

cm. (17) 

The proto-PN Hen 401 has cylindrical lobes with almost constant radius of 2;^ ~ 5 x 
10^^ cm up to a distance of ~ 3.5 x 10^^ cm from the central source (Sahai, Bujarrabal, & 
Zijlstra 1999). Farther out the lobes become narrower. Sahai et al (1999) present detailed 
study of the lobes. They argue, based on the theoretical study of molecular outflows by 
Masson & Chernin (1993), that the structure of Hen 401 can be explained by interaction 
between a well coUimated outflow and the slow wind, which intermediate between being 



~ 1 3 X 10^^ ( V^' ( X" 



tb{x) 
100 yr 



X 



obs 



1000 yr 



0.40 



X 



5 X W cm / V400 km s 



( 

V400 km s-V 
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radiative and being adiabatic. Some basic ingredients of their scenario apply here as well. 
According to the model presented here, the initial jet interaction with the slow wind is 
radiative, i.e., the shocked slow wind material has time to cool radiatively. This is appropriate 
for a jet blown by a main sequence companion; the case of a white dwarf companion is 
discussed below. As the jet expands to r ~ 10^^ cm (eq. 14 above) the adiabatic phase 
starts, with the inflation of a bubble (which will become the observed lobes; one lobe at 
each side of the equatorial plane) . This part of the bubble compresses the material near the 
equatorial plane; this explains the wide opening of the lobe close to the equatorial plane, 
which cannot be explained with a fully radiative flow. The jets then expand with a constant 
cross section. Further out, the jets may become weaker, and the interaction becomes unstable 
(Sahai et al. 1999). To quantitatively match the structme of Hen 401, 2:5 = 5 x 10^^ cm up 
to a; = 3.5 X 10^^ cm, the age of the jet and its velocity should be known. This information 
is not available, hence, 1 consider a few examples. We can take a^- = 3 x 10^^ cm, keeping the 
other parameters as in equation (17); or we can take tobs = 1500 yr, Vj = 500 km s^"*^, and 
keeping aj as in equation (17). For a white dwarf companion, we set Vj ~ 5000 km s^^ and 
take ~ 2 X 10^^ cm ~ 13 AU, keeping the observing time as in equation (17). In this case, 
the bubble inflation starts much closer to the binary system (eq. 14). To summarize, the 
flow structure presented here, despite its several assumptions, e.g., neglecting the decrease 
in the slow wind density along the z direction (perpendicular to the symmetry axis), can 
explain the structure of the lobes in Hen 401, and similar PNe having cylindrical lobes (see 
Sahai et al. 1999 for more examples). 



5. SUMMARY 

In the present paper, I studied the interaction of jets, or a coUimated fast wind (CFW), 
with a previously blown AGB wind, and reproduced some basic properties of bipolar FN 
lobes. In reality, the jets (or CFW) are blown by an accreting compact companion simulta- 
neously with the slow wind. To facilitate analytical solution, I had to assume axisymmetric 
flow, as well as several other simplifying assumptions. In addition, I studied only two ex- 
treme cases, a slow-propagating and fast-propagating jets, which occur when the slow wind 
density is much larger, or much lower, than the jet density, respectively. The goal was to de- 
rive simple expressions from which the properties of FN lobes blown by jets can be inferred. 
This can be used both to analyze observations (see discussions following eqs. 6 and 17), and 
to guide future numerical simulations. 

The main results can be summarized as follows: 
(1) The interaction of jets (or CFW) blown by main sequence stars, having a speed of 
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Vj ~ 300 — 500 km s~^, with the slow AGB stellar wind close to the central binary system is 
radiative, i.e., the shocked gas cools on a short time. Further away, the interaction becomes 
adiabatic, i.e., a bubble of hot gas which cools slowly is formed. The distance from the center 
where this transition takes place for the two types of jets studied here is given by equations 
(6) and (14). The expressions are scaled for jets blown by main sequence companions. Jets 
blown by white dwarf companions move much faster, Vj ~ 5000 km s~^, and start to inflate 
a bubble very close to the binary system. 

(2) As noted by SROO, the interaction between the CFW and the slow wind pushes material 
toward the equatorial plane. The formation of the inflated bubble at a distance from the 
equatorial plane in jets blown by main sequence stars, increases the efficiency of this process. 

(3) The significance of the transition from the radiative to the adiabatic phase to the struc- 
ture of the lobes, implies that radiative cooling must be incorporated in numerical codes 
intended to study the formation of bipolar PNe. 

(4) The basic interaction of a radially expanding jet with the slow wind forms an axisym- 
metric bubble which is narrow close to the central star(s) and far from the central star(s); 
its widest diameter is somewhere in-between. This type of structure is observed in many 
bipolar PNe. 

(5) If a jet is recoUimated, e.g., by magnetic fields inside the jet, such that its cross section 
stays constant, then the lobe widens close to the center, but then acquires a cylindrical shape 
with a constant diameter. Such lobes are observed in several bipolar PNe. 

I thank Liz Blanton for comments on the original manuscript. This research was sup- 
ported in part by grants from the US-Israel Binational Science Foundation. 
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